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ABSTRACT 



It has recently been argued that there are no dynamic screening 
corrections to Salpeter's enhancement factor in thermonuclear reactions, 
in the weak-screening limit. The arguments used were: 1) The Gibbs 
probability distribution is factorable into two parts, one of which, 
exp{—PJ2^i^j/Uj) W — ^/^bT), is independent of velocity space; and 
2) The enhancement factor is w = 1 + /3^e^ZiZ2(0^) with (02)^. = {E^)k/k'^ and 
(£;2)^/(87r) = (T/2)[l - £-1(0, A;)]. We show that both of these arguments are 
incorrect. 

Subject headings: nuclear reactions, nucleosynthesis, abundances 
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1. Introduction 

Our knowledge of reaction rates in the Sun is becoming more and more 
accurate. There have been several works dealing with electrostatic screening 
effects in the solar plasma (as for example, Carraro, Schafer, & Koonin 1988; 
Gruzinov & Bahcall 1997; Briiggen & Gough 1997). Recently, more precise 
calculations have been made beyond the linear regime (Gruzinov & Bahcall 
1998). 

The screening has the effect of lowering the Coulomb barrier between the 
interacting ions, therefore enhancing the reaction rates. This is included in the 
enhancing factor in the weak screening limit (where ZxZ2e^ j Rd^ 1 and 
is the Debye radius) 

w = exp A , (1) 

where 

A = Z^Z^e^/RoT . (2) 

This screeening factor uses the Debye-Hiickel expression. 

Usually, the calculation is made in the electrostatic case, based on the 
calculation of Salpctcr (1954). It is assumed that the motion of the screened ion 
is slow, compared to the motion of the screening particles. Carraro, Schafer, & 
Koonin (1988) have studied the case of dynamic screening. 

In a recent study, Gruzinov (1998) argued that in the weak screening limit, 
there is no dynamic screening corrections to Salpeter's enhancement factor, 
even for high energies. He based his conclusion on two arguments: 1) For 
thermodynamic equilibrium, the Gibbs probability distributions in velocity and 
configuration space arc decoupled; and 2) Through an analysis of the thermal 
electric field, we can estimate the random electrostatic potential and conclude 
that the enhancement of the reaction rates is given by the Salpeter expression. 
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We show below that although these arguments are simple, a careful 
examination shows that they are wrong. In §2 and §3 we recall the arguments 
and show why they are incorrect. 



2. Gibbs Distribution 

The enhancement factor for a reaction rate between nuclei of charges Z\e 
and Z-ie is 

w = exp{-Z2e<i)o/T) (3) 

{kB — 1), where 0o is the electrostatic potential, created by the plasma on the 
particle Zie. 

A test particle Zie moving through a plasma with velocity v' suffers dynamic 
screening. The electrostatic potential is written as (see Krall & Trivelpiece 
(1973), chap.ll) 



00 = 47reZi j 



(2vr) 



e(/c, kv') 



k-' , (4) 



where s is the dielectric permittivity, which describes the plasma response to 
the test particle, {e depends on the velocity distribution of the plasma particles 
f{v), which usually is taken as Maxwellian). Prom the above expression, it can 
be seen that 0o depends on the velocity v' of particle Z^e. Gruzinov (1998) 
argued that the Gibbs probability distribution p is 

(where /3 = 1/T) which can be factorable into 

p^exp(^-PY.^)^^p[-^i:^-^) ■ (6) 

From the above, it was argued that the distributions in velocity and configuration 
space are decoupled. 
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This first argument is simple. However, it is based on a misunderstanding. 
The general Gibbs probability distribution p for a plasma is 



where Wj^ is the interaction energy between the particles. Wik is the interaction 
energy of particle i with all the other particles in the plasma. This energy is the 
Coulomb energy eiCj/rij, related to the positions of all the other particles. It 
is an assumption that their positions Vj are independent of their velocities. 
However, we know from dynamic screening, Eq. (4), that their coordinates are, 
in fact, velocity dependent. Eq. (6) states that the coordinates are independent 
of the velocities only in zero order. The exact Gibbs distribution takes into 
account dynamic corrections. 

In fact, this argument is circular. The separability of the Gibbs distribution 
in velocity and configuration contributions is true only if the particles are 
moving under a (static) conservative force. The argument to show that there are 
no dynamical contributions, therefore, rests on a statement that is valid only if 
there are no dynamical contributions, which, in fact, is what he set out to prove. 

3. Thermal Electric Field 

Due to thermal fluctuations in the plasma, the reaction rate between two 
fast moving ions Z^e and is enhanced. The enhancement factor of a reaction 
rate is 



where (0^) is the average of the square of the random electrostatic potential 0. 
(0^) is given by 






(8) 




(9) 
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where (0^), = {E')Jk\ 

In his second argument, the expression used for the fluctuation electric fleld 
was (KraU & Trivelpice (1973), chap. 11) 



(Stt) 



1 - 



£(0,A;) 



(10) 



Substituing this expression in Eq.(8) and Eq.(9), Salpeter's expression is 
obtained. 

However, the expression of the thermal electric fleld (for a Maxwellian 
plasma), given by Eq. (10), assumes that a; <S T (^, fcs = 1). The 
general expression for the intensity of the electric field is given by the 
Fluctuation-Dissipation Theorem (see for example Sitenko (1967), Akhiezer et 
al. (1975)): 

(Stt) e-/^-l|£|2 ^ ' 

and 

^^'^^ - U.^'-^, . (12) 



(Stt) 



e-/^ - 1 I £ I' ■ 

This expression includes electrostatic Langmuir waves as well as all other 
fluctuations that exist in a plasma. In the limit of a; <^ T, 



(Stt) 



T Ime 



e 



|2 ' 



(13) 



for which it is then possible to use the Kramers- Kronig relations. Eq.(13) then 
turns out to be 

<^ = (L\ [l ^1 (14) 

which is the expression used. 

The assumption that uj <^ T, however, is very strong. In fact, in the case of 
the transverse electric fleld, we showed (Opher & Opher (1997a, 1997b)) that 
only by not making this strong assumption, is the blackbody at high frequencies 
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obtained. We recently showed (Opher & Opher (1999)), that by not making the 
assumption that a; <S T, the energy of a plasma in the classical limit is larger 
than previously thought. 

Without assuming that u <^T, the enhancement factor is given by Eq. (8) 
with Eq. (12). 

It is also to be noted that the second argument is also circular. It assumes 
that a; -C T, making it a static analysis, which is then used to prove that there 
does not exist a dynamic contribution. 

It is to be emphasized that we are not arguing here whether or not dynamic 
screening exists in thermonuclear reactions. We only show that the arguments 
used by Gruzinov (1998), to prove that dynamic screening does not exist, are 
not valid. 
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authors would also like to thank the Brazilian project Pronex/FINEP (no. 
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